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SYMBOLS 


A integral quantity (eq. (A3)) 

c speed of sound 

e total energy 
F radial flux 
/ dependent variable 

G axial flux 

// step height, ; or source of radial momentum (eq. (Al)) 

L location of pressure sensor 

M Mach number 

p pressure 

[/? ] Reynolds number in a cell, 

Rq base radius 
R^ sting radius 

R^ radius of the shock-tube wall 
r radial direction 

s distance from the root of the base: negative values indicate distance toward the shoulder; 
positive values indicate distances toward the downstream direction 

T temperature 

t dimensionless time, tHR^jc^ ) 

U free-stream velocity 

u velocity in x direction 

V velocity in r direction 

X axial distance from the root of the base 


hi 



7 specific heat ratio 

A dimension of a cell (size of cell, total length/ 127) 

5 increment 

Sfj Kronecker delta: 1 if / = / and 0 if / =5^ j 
e eddy viscosity 

f ratio of the reattachment distnace for 7 = 1 .4 to that for other values of 7 
p density 
j time 

CO ■ vorticity, co = - 
CO ^KciIRq) 

Subscripts: 

D dividing stream line 
i,i,k indexes 
o base 
R reattachment 
W wall 

00 undisturbed free stream 

1 behind the primary shock wave 


2 behind the bow shock wave 



STEADY AND NONSTEADY SUPERSONIC TURBULENT AFTERBODY FLOW 


Kenneth K. Yoshikawa and Alan A. Wray 
Ames Research Center 


SUMMARY 


A theoretical study of a supersonic turbulent afterbody flow field, containing regions of sepa- 
ration, recompression, recirculation, redevelopment, and inviscid flow, is made for axisymmetric 
supersonic steady and nonsteady flow past a cylindrical base supported by a sting. 

The effect of the specific heat ratio 7 of the incoming ideal gas on the flow properties, espe- 
cially on pressure distributions along the base and sting surfaces and on reattachment distance, is 
investigated. The specific heat ratios considered are 7 = 1.2, 1.4, and 1.667. Also, effects of other 
major parameters, such as eddy-viscosity coefficient (or effective Reynolds number) and Mach 
number, on the afterbody pressure and reattachment distance are studied and discussed. Evolution 
of shock-induced flow and stabilization time are examined and discussed for a transient problem. 

The important influence of the flow-field geometry, pressure distributions, and reattachment 
distance on the aerodynamics and radiative heat transfer for an atmosphere entry probe in high- 
speed flight are briefly described. 


INTRODUCTION 


The supersonic flow over a body with a blunt trailing edge (afterbody flow) is one of the most 
complicated problems in fluid mechanics. It contains a rapid expansion region around the base or 
trailing edge, a turbulent jet mixing of the external flow with a recirculating flow at the base region, 
recompression process, reattachment, and finally a redevelopment region due to the trailing shock 
wave at the wake. Because this afterbody flow strongly influences the drag of high-speed vehicles 
and missiles, considerable efforts have been devoted to its understanding (refs. 1-11). More 
recently, two other important projects involving the study of afterbody flow have been undertaken. 
These are the Galileo probe to Jupiter and the MX missile placed in a long ground tunnel. It was 
recently found that a significant amount of the radiation heating to the Galileo entry probe comes 
from the highly recompressed gas radiation from the wake region; this is in addition to the well- 
known radiant heating to the forebody surface (refs. 12-17). Thus, the design of the heat shield for 
the base of this vehicle requires a thorough understanding of the wake-flow structure from the near- 
wake to the far-wake regions. In the MX missile problem, the characteristics of the blast-induced 
shock wave and its interaction with the ground tunnel affect the survivability of the missile 
(ref. 18). This problem also involves the nonsteady turbulent flow over the missile, including the 
wall interference with the primary and secondary bow shock waves. 

Reviews of various methods of analysis applicable to the wake-flow regions have been pre- 
sented in references 10 and 11. Analytical methods using integral approaches (refs. 1, 2) and 



empirical models based on detailed experimental data (ref. 7) have been proposed previously. 
Recently, mathematical turbulence models have been reviewed in detail in reference 19. The turbu- 
lence models, however, are still not well understood. A simple model (of constant viscosity, con- 
stant Prandtl number, Pr = \, and constant 7) is desirable to simplify the computations; one such 
model has been previously introduced to solve the laminar near-wake-flow problem with some 
success (refs. 4, 5). 

Major emphasis in this report is on the development of a simple analytical method for studying 
the turbulent afterbody flow problem. A simple model, using an appropriate value of eddy viscosity 
and ideal gas, is employed for the fuUy developed turbulent flow calculations. This simple model 
will be particularly useful in providing some physical insight into the problem, for example, by 
examining the numerical effect of eddy viscosity (or effective Reynolds number) and specific heat 
ratio 7 on the flow field. 


NUMERICAL FLOW SIMULATION 


Steady Flow 

A schematic of the afterbody flow is shown in figure 1. The afterbody flow field consists of 
several distinct regions and processes: The supersonic flow, assumed uniform at the comer (with a 
thin boundary layer) at Mach number M-^ , expands rapidly and separates around the edge of the 

base, forms a free shear layer, and goes 
through a constant-pressure jet mixing, 
recompression, and reattachment; it then 
goes through a compression shock wave 
and redevelopment process. A portion of 
the shear-layer fluid is trapped in the base 
region and creates a recirculating flow 
bounded by a dividing streamUne, the 
base wall, and sting surface; the rest of 
the viscous flow proceeds in far-wake 
direction. 

To evaluate steady-state flow, two 
computer codes, the Weng-Chow code 
(ref. 1) and Wray’s code (ref. 20), were 
used, and some of their significant results 
are compared. Note that the Weng-Chow 
code has a more realistic turbulence 
model, but its solution is limited to flow 
in the recirculation region; on the other hand, Wray’s model is primitive (though it can be modi- 
fied), but its solution covers the entire flow region. Of the various fluid dynamics properties (pres- 
sure, density, temperature, velocity, vorticity, etc.), the pressure distribution is especially informa- 
tive at hypersonic speeds since it is mainly dependent on velocity and less so on the other fluid 
properties and physical parameters. These codes are briefly described below. 


SEPARATION EXPANSION RECOMPRESSION 



Figure 1 .— Separation and recompression regions of axi' 
symmetric flow over a cylindrical step base. 
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In the Weng-Chow code, the methods of analysis and calculation for various flow properties in 
the comer, turbulent jet mixing, recompression, reattachment, and redevelopment regions are 
obtained by the integral method described in detail in references 1 and 2. The inviscid flow region is 
described using the method of characteristics. The Prandtl mixing-length theory is used to compute 
the turbulent shear stress in the free mixing layer under a constant-pressure assumption. A spatial 
variation of eddy viscosity — in which the eddy viscosity is proportional to the value at the end of 
the mixing region times the local velocity scale and to the seventh power of the distance from the 
mixing region — was employed in the region of recompression. Solutions are obtained by cut-and- 
trial, adjusting some of the major parameters between trials. (The program was loaned to Ames 
Research Center by Professor W. L. Chow, University of Illinois.) The calculation of wall-pressure 
distribution along the sting surface was reported to show good agreement with experimental data, 
except that the computed peak pressures were generally higher and the gradients much steeper than 
experimental values (refs. 1 , 2). 

Wray’s code is a computational code for afterbody flow, using the compressible Reynolds- 
averaged, Navier-Stokes equations; it is an explicit method that is second-order accurate in space 
and time (unsplit, Richtmyer and Morton method). A brief derivation and description of the 
method are given in the appendix. A steady-state condition is reached in about 10,000 steps; it 
requires about 20 min on the ILLIAC IV. The computational domain considered is about 4 times 
the body radius in the radial direction and about 6.5 times the body radius in the streamwise direc- 
tion. Emphasis is on constructing the overall flow field around the afterbody. An adiabatic base and 
sting (and slip-flow condition for the shock-tube wall when wall interaction is considered) are used. 
A simple physical model has been incorporated in the calculations to determine the qualitative 
effects of the parameters on the flow-field 
solutions; the model consists of an ideal gas 
(constant 7), constant eddy viscosity, and 
unit Prandtl number {Pr =1). 

The essential features and results of the 
present calculations, using Wray’s code, are 
described in the following sections and are 
shown in figures 2-6. 


Effect of cell Reynolds number— The 
effect of cell Reynolds number [R ] , or 
inverse dimensionless eddy viscosity, on the 
flow field is shown in figures 2 and 3. Shown 
in figure 2 are computed pressure contours 
and streamlines in the afterbody region. 
Conditions for this case are [R] = 100, 
7= \2,Mi =2.38, andR^/R^ = 1/2. 

Conical expansion and recompression 
waves can be clearly identified in fig- 
ure 2(a). There are large constant-pressure 
regions behind the cylindrical base and 
downstream of the recompression shock 
wave. The slight variation of the pressure 



(a) Pressure contours. 



(b) Stream lines, exponential spacing. 

Figure 2.— Afterbody flow solution: Afj = 2.38, 
Rj/Ro = 1/2, [R] = 100, and 7 = 1 .2. 


3 



pd/pi 



x/Ro 


(a) Pressure contours. 



x/Ro 


(b) Stream function, exp (\p). 

Figure 3.— Afterbody flow solution: Mi = 2.38, 
1/2, [i?] =25, and 7 = 1.2. 



(a) The effect of [/?] on pressure distribution along divid- 
ing stream line: Mi = 2.09,7= l.67,Rs/Ro= 1/2. 

Figure 4.— Pressure distributions. 


between the base and the recompression 
shock indicates a possibility of boundary- 
layer-like flow in that region. The pressure 
distribution along the sting surface shows 
that the pressure buildup is gradual; after the 
peak value is reached, the pressure decreases 
very slowly to the free-stream value. Shown 
in figure 2(b) are streamlines (with exponen- 
tial spacing); a pattern of the recirculatory 
flow as well as the dividing streamline 
through the rear stagnation point (reattach- 
ment points are marked by circular symbols) 
can be seen. 

Pressure contours and streamhnes for a 
larger eddy viscosity ([/?] = 25) are shown 
in figure 3 for the same values of y, Mi, and 
R^/Rq as used in the previous case. The base 
pressure below the comer is drastically 
reduced (fig. 3(a)). Strong vorticity and tur- 
bulence carry momentum around the comer 
and delay the boundary-layer separation, 
that is, separation occurs below the comer, 
decreasing the size of the wake and the 
reattachment distance, as can be seen from 
figure 3(b). This delay of boundary-layer 
separation has been observed in numerous 
experimental investigations (e.g., ref. 9). 
Some numerical analyses for laminar two- 
dimensional flow behind a wedge or step 
(refs. 4, 5) also pointed out this behavior. 
The separation delay will be discussed fur- 
ther in a following subsection (Vorticity 
Distribution). The steep pressure gradient 
near the comer suggests that the boundary- 
layer character of the flow along the comer 
vanishes as the flow expands. These results 
simply indicate the influence of larger eddy 
viscosity. In our report, however, for sim- 
phcity we assume a constant eddy viscosity. 

The effect of [/? ] on pressure distribu- 
tions is shown in figures 4-6. Pressure distri- 
butions on the dividing streamline, particu- 
larly in the vicinity of the comer, vary 
significantly with the value of [/?] , as shown 
in figure 4(a). The results for other values of 
7 show pressure variations similar to those 
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(b) Pressure distributions for the effective [/?] =75. 
Ml =2.09,7= 1.67, 1/2. 


for 7 = 1 .67. Note that the pressure distribu- 
tion along the sting surface does not 

vary as much with [/?] as that along the 
dividing streamline (p£))- Also note that a 
steep pressure gradient appears near the 
comer for the small value of [i?] ; this is also 
clearly shown in figure 3(a). 

The concept of constant-pressure mix- 
ing, after the flow separation, has been sup- 
ported by many authors for lower super- 
sonic speeds, that is, <3 (e.g., refs. 1-6, 

10,11, 13). In order to obtain relatively 
constant pressure in the base mixing 
region and also to keep the base pressure 
close to this pressure (i.e., pg = pj^), an 
appropriate value of [R] has to be selected. 

The value that satisfies these conditions is 
[R] ^75. Pressure distributions of pj^ and 
p^; for [7^1 =75 are shown in figure 4(b) 
where, as can be seen, P/) in the turbulent 
jet mixing layer is nearly constant from 
x/fl = 0 to xjH 2.0, ahead of the reattach- 
ment point, and is also nearly equal to the 
base pressure (note that p^ becomes identi- 
cal to after flow reattachment). 

The numerical calculations for other 
values of 7 (i.e., 7 = 1 .2 and 1 .4) show simi- 
lar results and, thus, support the choice of 
[R] =75. Presented in figure 4(c) are pres- 
sure distributions for a higher Mach number 
(Ml = 4). The numerical results show that 
the effect of [R] on pressure distributions is 
relatively small for higher Mach numbers, 
but is somewhat larger for lower Mach num- 
bers. Thus, base pressure is weakly depen- 
dent on [R] (see also refs. 1, 10). It is also noted that pressure varies rapidly along the dividing 
streamline and that constant-pressure mixing can no longer exist in this high-shear flow region for 
higher Mach number. Note that the feature of constant wall pressure (below the separation point) 
still exists in the base flow region. 



(c) The effect of [R] on pressure distributions at higher 
Mach number. Ml = 4.0 , 7 = IA,R^/Rq = 1/2. 

Figure 4.- Concluded. 


In figure 4(c) negative values on the abscissa indicate locations between the comer (s/H = -1) 
and the root of the base (s/H = 0 ). 


Effect of specific heat ratio— For hypersonic entry into a planetary atmosphere it is known 
that the specific heat ratio of the high-temperature real gas will change significantly during the 
flight; therefore, the examination of the effect of 7 on the flow field is important. Three different 
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(a)Mi = 2.09, RsiRo = 1/2, and [^] = 75. 



f • s/H 


values of constant y are used in the calcu- 
lations in order to see the effect of y on 
the pressure distributions. 

A rather good correlation of pressure 
distributions for 7 = 1.2, 1.4, and 1.67 is 
obtained in the form of a reduced coordi- 
nate, as shown in figure 5(a). The flow con- 
ditions for the cases are Mi = 2.09 and 
Rg/Ro ~ 1/2. The factor f in the abscissa is 
the ratio of the reattachment distance for 
7 = 1 .4 to that for a given 7. Experimental 
correlation of pressure distributions in the 
reduced form have been demonstrated in 
reference 6 for air, and seem valid for a rea- 
sonable range of R^/Ro- The detailed effect 
of the ratio of the sting radius to the body 
radius, RJR^, on the pressure distributions 
and reattachment distance are shown in 
references 1 and 8, based on theory and 
experiments in air. It should be noted here 
that significant overshoot pressure, following 
the flow reattachment, is observed for axi- 
symmetric flow (see also fig. 4), but no over- 
shoot is seen in two-dimensional flow 
(refs. 6, 8). 


(b)jWi -4.0, [/?] - 75, and^j//?o - 1/2. jhe effect of 7 on the pressure distri- 

Figure 5.- The effect of 7 on pressure distributions. bution at a higher Mach number, M 1 = 4.0, 

is shown in figure 5(b). The results clearly 
indicate that the base pressure is strongly dependent on both flow Mach number Mi and specific 
heat ratio of the gas 7, that is, the base pressure ratio is significantly reduced by increasing and 
by decreasing 7. The maximum pressure ratio, Pmax/Pi > however, is relatively independent of these 
parameters (Pmax/i^i slightly increases asA/j increases). Consequently, pressure distributions in the 
reduced coordinate do not show as good a correlation as for the lower Mach number case in the 
region of near-wake, but retain good correlation in the far-wake region. In figure 5 negative abscissa 
values indicate the locations from the comer (s/// = -1) to the root of the base {s/H = 0). 


Vorticity distributions— Vorticity is plotted in reduced coordinate in figure 6 for the higher 
Mach number condition. The vorticity along the dividing stream line (oj£)) and along the base-to- 
sting surface (copj/) is presented. Vorticity is large around the comer. Three zero vorticity points 
(cj = 0), one just below the corner (an evidence of the delayed separation point), one at the base 
(rear stagnation point), and one at the flow reattachment location are present. Note that the vortic- 
ity oi£, in figure 6 should start from the origin, that is, = 0 at s = 0. The vorticity distributions, 
W£) and ojpj;, become identical after flow reattachment. In contrast to the pressure distributions in 
figure 5(b), good correlations of the vorticity distributions exist in the region of recirculation of 
near-wake flow. The vorticity along the sting and base surface is positive inside the recirculation 
region, corresponding to reverse flow; along the dividing stream line outside the recirculation region 
(along the surface before the separation and after the reattachment) it is negative, corresponding to 
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flow in a streamwise direction. Note that toyj, along the sting surface 

strong vorticity is regenerated as the flow g “d along the dividing streamline 

moves past the reattachment point. 

"w 

Other general physical results are ^ 

described in the following section — General 

Results and Discussions. 4" 
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Nonsteady Flow 

A shock tube is a useful facility for the 
experimental study of high-speed entry gas ? * 

dynamics (ref. 21). In particular, shock-wave 
interaction with walls and projectiles inside a 

shock tube or in a long ground tunnel is an Figure 6.- The effect of y on vorticity distributions: 

important subject of investigation for heat =4.0, [R] =75,and/?j/RQ = 1/2. 

protection and safeguarding of space vehicles 

and missiles (ref. 18). Numerous reports have presented various features of numerical calculations 
for the inviscid shock-induced flow, that is, supersonic wake flow (ref. 22), shock-wave interaction 
with the cross-flow (ref. 23), and flow patterns of diffracting weak shock waves around comers 
(ref. 24). The important nonsteady flow characteristics associated with a shock-tube test, such as 
actual test time, flow stability time, and flow reflection time, are not well understood, except that 
the reflection time (time for the shock wave to travel to the end wall and return) is usually much 
longer than the test time. In this section the effects of turbulent flow on the transient pressure, flow 
stability time, and some effects of 7 on the pressure distribution in the shock-tube flow will be 
briefly described. 

Shock-wave-induced base flow— A sketch of the initial shock wave passing over a sting- 
supported cylindrical base is shown in figure 7. The primary shock wave at Mach number Mqo dif- 
fracts at the corner and passes to the sting 
surface, developing the turbulent shear 
layers. The principal results for the transient 
phase of flow, from nonsteady to steady- 
state conditions, are shown in figure 8 for 
the condition of = 6, 7 = 1.4, and 
Rs/^o “1/3- Figure 8(a) shows the density 
contours and pressure distribution along the 
sting surface shortly after collision of the 
diffracted shock wave with the base line 
(f = 0.37) and after collision with the sting 
surface (F= 0.74). The pressure distribution 
at 7 = 0.37 is still undisturbed, since the 
shock-wave front has not reached the bot- 
tom of the base surface (or sting surface). 

Strong overshoot pressure, due to shock 

wave reflection or diffraction, develops Figure 7.— Geometry of shock wave passing over a 

immediately and reaches a maximum value cylindrical step base supported by sting. 
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(a) Density contours at / = 0.37; (b) Density contours at / = 1.49; 

wall pressure at 7 = 0.74 and 0.37. waU pressure at 7= 1 .09 and 1 .49. 

Figure 8.— Evolution of density and pressure distributions for shock -induced base flow: 

= = 1/3, [7?] =75,and7= 1.4. 

of pip I « 2.9 (where pj is the pressure behind the incident shock) at 7= 1.0, as shown in the subse- 
quent figures. The major portions of flow field follow the inviscid flow character following the 
conical expansion around the comer, until flow time reaches t~ 1.0; no significant damping or fluc- 
tuation due to viscous effects appears in the pressure distributions (see e.g., refs. 22, 23). Strong 
vorticity is generated as the shock moves past the comer. The starting vortex moves downstream, 
and its interaction with the sting becomes important, as can be seen in figure 6. 

The maximum pressure buildup behind the shock occurs at a characteristic time at F = 1.0 
without regard to the gas property (or 7) used in the numerical calculations. Radiant luminosity 
believed to be due to this peak pressure has been observed in reference 15. 

Flow stability time— The base pressure ratio (p/pi at xlR^ = 0.0) also varies with time; this 
effect can be seen in figure 8. Time history of the base pressure in terms of the reduced time is 
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shown in figure 9. The first pressure peak 
(impact pressure) occurs immediately after 
the diffracted shock wave impinges upon the 
sting surface at the bottom of the base, and 
the base pressure then builds up again, with 
some fluctuations, to the asymptotic value. 
We arbitrarily define flow stability time t 
by 

P5W/P'5(°°) = 0-90 (1) 

From several numerical calculations, the fol- 
lowing empirical relation is established for 
the case of =1/3, 

r = 61-(V^i)/(^i/'"oc) (2) 

The flow stability time will be dependent on 
the ratio of the sting radius to the body 
radius, R^/R^ \ it will be shorter when R^IRq 
is smaller since interaction with the bound- 
ary layer is less. 

It is possible that the flow stability 
time given by equation ( 2 ) can be longer 
than the actual flow test time, since the 
flow test time can be significantly shortened 
by the shock-tube boundary layer growth 
and by the broadening of the mixing regions 
across the contact surface in hypersonic 
shock tube. It is possible that the steady- 
state flow may not be established in many 
test cases. 



0 2 4 6 


x/Ro 


(c) Density contours at r = 115; wall pressure at 
r = 2.88 and 115. 

Figure 8.— Concluded. 


Shock-tube flow simulations— Typical results of the shock-tube flow simulation (F » 1) are 
shown in figure 10. The pressure distributions along the sting surface are plotted in this figure for 
three different gases (i.e., 7 = 1.2, 1.4, and 1.67) - 14.0 and R^IRq “ differ- 

ences in base pressure ratios are due to the different values of the specific heat ratio 7 , which also 
result in different flow Mach numbers Mi (Mi = 2.80, 1.86, and 1.33 for 7 = 1.2, 1.4, and 1.67, 
respectively); however, the peak pressure ratio Pmax/Pi not change significantly. The rapid 
pressure decay shown for a monatomic gas ( 7 = 1 .67) is due to shock-tube wall interaction (reflec- 
tion of the expansion wave from the shock-tube wall), a result of the lower supersonic speed of 
Ml = 1.33. The wall-interaction effect on pressure for other gases (7 = 1.2 and 1.4) would appear 
further downstream; it cannot be seen in these figures because the Mach numbers are too high. 
Validity of the concept of constant-pressure mixing in the recirculation region (below the separa- 
tion point) for various 7 and Mach numbers Mi is well substantiated in the figure, that is, the 
constant-pressure mixing and constant base-pressure assumptions seem to be valid for 7 > 1 .2 or for 
Ml <3.0. 
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Figure 9.— Transient base pressure and flow stability Figure 10 — Pressure distributions and wall interference 
time: boundary conditions same as in figure 8. (shock-tube flow simulation): Moo= 14, Rg/Rg —1/3, 

and Ry,/R() = 4. 

GENERAL RESULTS AND DISCUSSIONS 


Important physical results, such as the ratio of the neck pressure (maximum pressure) to the 
base pressure, P^iax/Pi ’ reattachment distance xj^, are presented and discussed in this section. 
(Note thatpjj^ax/Pi = (Pmax/P^) * (Pb^pO-^ 

Comparison of the Ratio of Base Pressure to Neck Pressure 

This pressure ratio is important in the estimation of the magnitude of radiation heat from the 
wake (see, e.g., refs. 15, 16). The computed peak-to-base (or neck-to-base) pressure ratios are shown 
in figure 1 1(a). Since the range of sting-to-body radius ratio in this report is limited to Rg/Rg = 0.33 
and 0.5, the radius effect on base pressure is relatively minor. Significant effect of flow Mach num- 
ber and 7 on Pf^ax^P^ noted as Mi exceeds 3.0. One can, therefore, expect very high pressure 
ratios for the hypersonic flow where dissociation and ionization of the gas begins to take place, 
that is, for a gas of low 7 . Also shown in the same figure, with solid black symbols, are the results 
obtained by using the code of reference 1 . An empirical relation derived for 7 = 1 .4 is presented in 
the figure as a dotted line, which agrees well with the results of reference 1 . 
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(a) Ratio of peak pressure to base pressure 
as a function of Mach number. 


(b) Base pressure as a function of 
Mach number. 


Figure 11.— Wake pressures. 


Base Pressure 

Base pressure ratios, p^lpi, calculated by the present code and by that of reference 1, are 
shown in figure 11(b) (numerous experimental and theoretical data are also presented in ref. 7). 
Significant effects of y on the ratio is noticed as Mach number Mi increases. The large effect of y 
on the peak pressure ratio shown in figure 1 1 (a) is due mainly to the effect of y on base pressure 
since Pniax/^i relatively unchanged (see figs. 4, 5). The result shows that base pressure will be 
significantly lower when y is lower, that is, when dissociation or ionization takes place in hyper- 
sonic flight. Note that the ratio of sting radius to body radius is also limited to RJRq = 0.33 
and 0.5 in the figure (the effect of body radius on base pressure, as mentioned before, is demon- 
strated in references 1 and 8 for lower supersonic Mach number). 


Reattachment Distance 

Reattachment distance of the afterbody, shown in figure 12, is an important characteristic 
length. It provides an approximate location of the vertex of recompression shock and, thus, gives an 
appropriate view factor for the wake radiation involved with the high-temperature gas. Also shown 
in the figure is an empirical relation for the reattachment distance for y = 1 .4. A few computed 
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Figure 12 — Flow reattachment distance as a function 
of Mach number. 

viscous flow over a blunt forebody at supersonic 
can be incorporated with the present afterbody c 


results from reference 1 are also included in 
figure 12, but the results show strong depen- 
dence on 7 and fluctuate with a wide varia- 
tion (the method, however, is not expected 
to be applicable for higher Mach numbers). 
The present results indicate rather good 
correlation, whereas those from reference 1 
do not. 


Forebody Flow Effect 

The presence of a forebody with 
accompanying bow shock wave may signifi- 
cantly affect the afterbody flow solution. 
For example, the trend of experimental 
pressure distributions downstream of the 
separation point for various forebody shapes 
(and sting radius to body radius ratio) has 
been shown in references 6 and 8: distinct 
overshoot pressure, as shown in previous 
figures, results for the long cylindrical base 
flow, and no significant overshoot pressure 
is reported for the particular forebody pre- 
sented. Note also that the result of two- 
dimensional flow (or Rg/Rg — 1 .0) does not 
show the overshoot pressure. The detailed 
study of the forebody effect on the after- 
body flow is, therefore, desirable. There are 
a number of time-dependent finite- 
difference Navier-Stokes codes that have 
been developed to calculate the turbulent 
transonic speeds (e.g., refs. 25, 26). These codes 
; in order to see the effect of a forebody. 


Code Limitation 

Besides the flow-model limitations mentioned earlier (i.e., ideal gas, constant turbulent eddy 
viscosity, etc.), there is another code limitation. Solution of the shock-induced flow (nonsteady 
shock-tube flow and steady wind-tunnel afterbody flow) involves a steep pressure jump across the 
shock wave (including the recompression shock wave). These discontinuities often cause numerical 
problems. The relatively gradual pressure rise across the recompression shock wave shown in fig- 
ure 4, and across the diffracted shock wave shown in figure 8, can be attributed to this code limita- 
tion. (However, higlier eddy viscosity also causes a gradual pressure rise.) For example, pressure dis- 
tributions calculated by the code of reference 1 show much steeper pressure rises than the present 
results at the beginning of the recompression region, and the experiments seem to support the 
results of reference 1 , except that the later portion of the experimental pressure rise is much more 
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gradual than that predicted by reference 1 . 
Therefore, further examination and improve- 
ment on this subject may be essential. 


Shock-Tube Experiments 

To demonstrate the general trend of 
unsteady afterbody flow, a typical result of 
a shock-tube experiment for transient pres- 
sure distribution is shown in figure 13. The 
test model is a hemispherical body sup- 
ported by a sting, with a ratio of sting radius 
to body radius of R^/Rq = 1/3; the test gas is 
air at a shock-wave Mach number of about 
Mao ~ 14. As noted in the previous section, 
the forebody shape alters afterbody pressure 
distributions, especially after the reattach- 
ment point, and for the demonstration pur- 
pose we choose transient pressure near the 
base {x/H = 1.0), which seems to be less 
affected. Air at this test condition is par- 
tially dissociated, and two numerical calcula- 
tions (one for 7 = 1.2 and the other for 
7 = 1.4) are shown for the comparison. After 
a sharp pressure spike, due to the diffracted 
shock wave compressing and interacting with 
the boundary layer near the surface, the 
transient pressure gradually increases toward 
the asymptotic values of Pi, the steady-state 
pressure behind the primary shock wave. In 
this test condition we estimated that the 
experiment test time (300-400 psec) and 
the flow stability time (t = 300 /usee for 
7 = 1.4, and t = 400 /usee for 7 = 1.2) were 
comparable. However, it could be that the 
flow does not reach the steady-state condi- 
tion before the mixing region of the inter- 
face arrives at the location of the pressure 
sensor. Note again that the numerical calcu- 
lation of transient pressure is based on the 
infinitely long forebody and thin-boundary- 
layer assumptions; agreement of the test 
data and numerical calculations is not 
emphasized here. 



(a) 7 = 1 . 2 . 


1.0 r- 
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Q. 





50 /Ltsec 
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u 1 1 1 1 

0 t 10 20 30 40 

’■ TIME,t/(Ro/ci) 


(b) 7 = 1 .4. 


The extremely strong and nonsteady 
pressure induced by the blast shock 


Figure 13.— Experimental and theoretical transient 
pressure: = 14, and = 1/3. 
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wave, which hammers the body and thus affects stability and safety of the vehicle, should be 
noted. 


CONCLUDING REMARKS 


A computer code and solutions of steady and nonsteady supersonic turbulent-wake flow have 
been obtained under simple flow model assumptions. The emphasis of the study was on examining 
the basic characteristics of the afterbody flow field to obtain essential flow information for the cal- 
culation of radiation heat from the wake region. Solution of the entire afterbody flow covered the 
base through far-wake regions; within the framework of the present assumptions, the solution 
results in the following general conclusions: 

1 . There is a steep pressure jump (thus, high temperature jump) across the recompression 
shock wave (which can be correlated for low Mach numbers, < 3). 

2. Relatively large constant-pressure zones exist. 

3. There is a gradual decay of the pressure distribution in the far-wake region (thus, slow 
decay of temperature and radiation intensity). 

4. A conical shape of high-pressure zones (thus, high-temperature radiation sources) is an 
appropriate geometry for calculating radiation heat from the wake region. 

For nonsteady flow, the following results are noted: 

1 . There is a very high induced pressure on the sting surface, which is also substantiated by 
experiments (several times higher than pi , steady-state pressure behind the incident shock wave). 

2. Delayed separation, rear-stagnation, and reattachment points can be identified by the vor- 
ticity distribution curve. 

3. Flow stability time may be comparable with flow test time. 

However, further examination of the wake flow — using improved models that include, for 
example, the effects of real gases, variable turbulent viscosity, presence of forebody flow, and 
improved resolution across shock waves — is essential to the accurate prediction of radiation heat 
and heat protection of a very high speed entry probe, such as that for the Galileo mission to Jupiter. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, California 94035, July 23, 1980 
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APPENDIX 


DESCRIPTION OF WRAY’S CODE 


The Wray code solves the Navier-Stokes equations in an axisymmetric geometry by means of a 
finite-difference algorithm of second-order accuracy in space and time for the inviscid terms and 
second-order accuracy in space and first-order accuracy in time for the viscous and heat conduction 
terms (ref. 20). The method of differencing enforces exact numerical conservation of mass, momen- 
tum, energy, and angular momentum. 


We can write the Navier-Stokes equations in the following shorthand notation: 

a/j. 1 3F. ac.- 


(Al) 


where = (p, e, pu, pv) is an array of the dependent variables: mass density, total energy density, 
radial momentum density, and axial momentum density. and are the corresponding radial and 
axial flux densities. The term H, which occurs only in the radial momentum equation, constitutes 
a source of radial momentum due to geometric effects in the axisymmetric coordinate system. The 
operator is the Kronecker delta = 1 if / = 3, zero otherwise). 

We take node points rj = (/ - 1) 6/- and = {k - 1) 8x,j, k = 1, 2, . . ., 128, to represent the 
flow field. The values of the /j at (ry, x^) are considered to be average values in a cell. We integrate 
(Al) to find these averages: 


i A 

A dt 



(A2) 
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where 



The integrals in equations (A2) and (A3) are approximated by the trapezoidal rule. This constitutes 
a so-called “finite volume” approximation. 


The at a time t are first advanced inviscidly to a time /+ 1/2 5/ at the cell centers irj+ 1/2 5r, 
+ 1/2 bx), and then the inviscid flux densities, as computed at those points, combined with the 
viscous fluxes as computed at time t, are used to advance to time t + bt. 


Boundary conditions used are as follows. At the incoming boundary, where the flow is super- 
sonic, all dependent variables are specified as desired. At the downstream (outflow) boundary, a 
zero gradient (df^ldx = 0) condition is used; this condition is also used on the top boundary (as 
bf^ldr = 0) in the case of no outer wall. If there is an outer wall, slip boundary conditions are used 
(v = 0); on the body surface (and sting surface, if present), a no-slip, adiabatic condition is used. 
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